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Abstract. 

We study the Renyi entanglement entropies of two disjoint intervals in XY 
chains. We exploit the exact solution of the model in terms of free Majorana 
fcrmions and we show how to construct the reduced density matrix in the spin 
variables by taking properly into account the Jordan-Wigner string between the 
two blocks. From this we can evaluate any Renyi entropy of finite integer order. 
We study in details critical XX and Ising chains and we show that the asymptotic 
results for large blocks agree with recent conformal field theory predictions if 
corrections to the scaling are included in the analysis correctly. We also report 
results in the gapped phase and after a quantum quench. 
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= T ^\ogTrp a A , (1) 



1. Introduction 

The bipartite entanglement for a given division of the Hilbert space into a part A and 
its complement B, can be measured in terms of the Renyi entropies [T] 

1 

1 

where pa — Tr^ p is the reduced density matrix of the subsystem A, and p = \^){^\ 
is the density matrix of the whole system in a pure state jvP). The knowledge of the 
for different a characterizes the full spectrum of non-zero eigenvalues of pa [2] , 
and gives more information about the entanglement than the more commonly studied 
von Neumann entropy S^K It also gives a fundamental insight into understanding 
the convergence and scaling of algorithms based on matrix product states [5] . 

In [H [51 [3] it has been shown that for a one-dimensional critical system whose 
scaling limit is described by a conformal field theory (CFT), in the case where A is 
an interval of length I in an infinite system, the asymptotic behavior of the quantities 
determining the Renyi entropies is 

Trp2~c Q ^-J , (2) 

where c is the central charge of the underlying CFT. Thus the Renyi entropies (and 
in particular the von Neumann one for a = 1) give one of the best way of detecting 
the value of the central charge. 

Less attention has been devoted until now to the entanglement of two disjoint 
intervals in a CFT (and also in massive theories). It turned out that entanglement of 
disjoint intervals is sensitive to universal details of the CFT that are not encoded in 
the central charge and it is connected with the full spectrum of operators of the CFT 
underlying the lattice model. 

We consider here the case of two disjoint intervals A — Ai U A 2 — [tti, v{\ U [u 2 , v 2 ] 
depicted in Fig. [T] By global conformal invariance, in the thermodynamic limit, Tr p°^ 
can be always written as 

/ I II I \ — 

rr, a 2 I Ml -U 2 \\Vl -V 2 \ V 771 / \ r „N 

Ti pa = c a I n n n , F a ( x ) , 3) 

VK - Vi\\u 2 - V2\\Ul - v 2 \\u 2 - Vi\) 
where x is the four-point ratio (for real Uj and Vj, x is real) 

x = (ui -vi)(u 2 -v 2 ) 
(ui - u 2 )(vi - v 2 ) ' 

Normalizing such that F a (0) — 1, we have that c a is the same non-universal constant 
appearing in Eq. The universal function F a (x) depends explicitly on the full 

operator content of the theory and must be calculated case by case. Originally it was 
proposed that F a (x) = 1 identically [5]. This erroneous prediction has been tested 
in free fermion theories [71 [51 IH| , and only the calculations in more complex theories 
allowed to detect this error [TU1 [TT1 [T2 l [TB ] . 

Furukawa, Pasquier, and Shiraishi [TT] calculated F 2 (x) for a free boson 
compactified on a circle of radius R 

P / x 9 3 {vt)9 3 (t/ v ) 
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where r is pure-imaginary, and is related to x via x = [02 (t) / '^(t)] 4 . 6>„ are Jacobi 
theta functions, rj is a universal critical exponent dependent on the compactification 
radius R (in Luttinger liquid literature r\ = 1/(2K)). This has been extended to 
general integer a > 2 in Ref. [H] 

e(o|»?r) e(o|r/r?) 

where T is an (a — 1) x (a — 1) matrix obtained in [12 
is the Riemann-Siegel theta function 

9(0|r) = exp [in m t ■ T ■ m] 

m G Z°- 



F a (x) 



(6) 

r\ is the same as above, while 



(7) 



The XX model considered in the following is described by the compactified boson 
with rj — 1/2. We mention that the analytic continuation of this result to real a for 
general values of rj and x (in order to obtain the entanglement entropy) is still an open 
problem, but the results for x <C 1 and rj <C 1 are analytically known [12"] . 

Nowadays, the only other example of non-trivial F a (x) known exactly is F2(x) 
for the Ising model that has the rather simple expression |13| 

1/2 



F 2 (x) 



1 



(l + V5)(l + yT^) 

2 



1/2 



+ x 1/4 + ((l-a;)a;) i/4 + (l-x) 



1/4 



a/4 



.(8) 



We mention that all these new results for entanglement entropy were derived by using 
some old 'classical' results about CFT on orbifolds fPf] . 




Figure 1. Typical bipartition we consider in this paper. The subset A is the 
union of two disjoint intervals Ai and A2 of length l\ and £2 respectively. The 
block separating them is denoted with B± of length r. The 'environment' is 
B = B\ U B2- The thermodynamic limit is obtained by sending the total length 
L — > 00, while ^1,^2)'" remain finite (i.e. the length of B2 goes to 00). 
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It is worth to recall that in the case of many intervals, the entanglement entropy 
measures only the entanglement of the disjoint intervals with the rest of the system. 
It is not a measure of the entanglement of one interval with respect another, that 
instead requires the introduction of more complicated quantities because A\ U Ai is 
in a mixed state (see e.g. Refs. [T]3 [13 [H] for a discussion of this and examples). 

In this manuscript we report the exact evaluation of the Renyi entropies for integer 
a in XY spin chains with periodic boundary conditions, whose Hamiltonian is given 

by 

N 1 1 _ , 

Hxy = -J2 [^of of +1 + ^afaf +1 + £«f] , (9) 
1=1 

where <r" are the Pauli matrices at the site I, h is the magnetic field and 7 is the so- 
called anisotropy parameter. For 7=1 the Hamiltonian reduces to the Ising model, 
while for 7 = to the XX model. Hamiltonian ^ can be diagonalized in terms of 
free-fermions (see below) but the fundamental difficulty is that the reduced density 
matrix of spin and fermion variables are not the same (T3J [IS] ■ While calculating the 
fcrmion matrix is straightforward, we derive here a general method to deal with spin 
reduced density matrices and present actual calculations for the Hamiltonian ([9]). 

The paper is organized as follows. In the next section [2] we give a general 
description of the density matrix mapping and we introduce our formalism. In section 
[3]we derive the expressions for the Renyi entropies of integer order. These two sections 
are rather technical and are intended only for those readers that want to repeat the 
calculations. The readers interested only in the results can skip directly to section [4] 
where we report our results for the XX chain and compare them with CFT results by 
properly detecting corrections to the scaling. In section[5]we present the same analysis 
for the Ising model and we check the universality while changing 7. In section [6] we 
consider non-critical chains. In section [7] we consider the Renyi entropies after a 
quantum quench. Finally in section [8] we summarize our main results and we discuss 
topics deserving further investigations. Four appendices contain technical details of 
the calculations and some background material. 



2. The reduced density matrix 

Generically a spin 1 /2 chain can be described by means of the Pauli matrices tr^ 1 with 
/Lt = 0,1,2, 3 ((7° = 1 and 1 = x, 2 = y, 3 = z) and i labels the lattice sites. We are 
interested in the entanglement of a subsystem A = U*=i Aj consisting of N disjoint 
spin blocks Aj = [uj, Vj]. The reduced density matrix (RDM) can be written as a sum 
over all operators belonging to the blocks forming A 

^=2wE<ik)ik- d°) 

Mi leA leA 

where \A\ = J^j 1S the sum of the lengths of the blocks £j. The expectation value is 
taken over the state we need to calculate the entanglement. The knowledge of multi- 
point correlation functions ^H;eA (7 f ' ) ^ s nowever a very hard task even for eigenstates 
of integrable models (see e.g. [TH1 HO] for the XXZ chains) and the above formula is 
particularly useful only when the spin-chain admits a representation in terms of free 
fermions, for which the Wick theorem suffices to calculate all of them. 
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The mapping between spin variables and fermion ones is obtained via the Jordan- 
Wigner transformation (that for the XY chains makes the Hamiltonians quadratic) 

-=(n<)^> 4=(n<)^- (id 

It is also useful to define the Majorana fermions [ST] 

a 2 i=cj + a, a,2i-i = i(ci - c\) , (12) 
satisfying anti-commutation relations {a/, a n } = 2Si n . For the case of a single interval, 
a fundamental observation 21] is that the Jordan- Wigner string frim<J a m) ma P s the 
space of the first £ spins in the space of the first £ fermions and thus the RDMs of the 
first £ spins and fermions are the same, i.e. 




PA = ^ > K' < , (13) 



with — 1 and aj — a;. At this point, it is worth to remind that the above equation 
gives directly all the eigenvalues of pA and so all Renyi entropies (see e.g. the reviews 
[22] for more details). Indeed, one recognizes pa as an exponential form 

p A oc e a ' w '- a -/\ (14) 



and using the Wick theorem (see also Appendix A ) 
W 

tanh — = r, with Tij = {cLidj) — Sij . (15) 

The 2£ eigenvalues of the correlation matrix T have the form ±z/j with vi £ [— 1, 1], so 
one has all the 2 e eigenvalues of pa as 

Am - II ■ ( 16 ) 

i 

for the 2 e possible choices of the variables Si — ±1. The Renyi entropies for any 
complex a are then given by 



1 + vA f l 



(17) 

J = l 

The calculation of vi only requires the diagonalization of a 2£ x 2£ matrix, that is 
an extremely easy task compared to the diagonalization of the full 2 l x 2 l density 
matrix. This expression can be further analytically massaged to obtain exactly the 
non-universal constant c a and also universal corrections to the scaling pM] . 

This is no longer true in the case of more disjoint spin-blocks. In the fermionic 
space, the subspace of two disjoint spin-blocks is not a local quantity: the non-locality 
of Jordan- Wigner transformation causes the fermions between the two blocks (i.e. 
those in Bi) to contribute to expectation values. For example, considering the product 
OjO-% with I 6 Ai and n £ A 2 and calling r the distance between the blocks (see Fig. 
FT]), we have 



(18) 
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that depends on the 'unpleasant' string of Majorana operators 

ti+r 

S= Yl ( ia 2j-ia 2 j) ■ (19) 

The string contribution cannot be neglected (as also shown in [T3JIIH]) and so the spin 
representation is not equivalent to the fermionic one. 

For practical reasons we limit ourselves to the case of two blocks \A%\ = 1% and 
\A%\ = 1%. We denote with \B\\ = r the space between the two blocks, see Fig. [I] for 
a graphical representation. We focus on spin-chains with the ground-state property 
that only even numbers of fermionic operators have non-vanishing mean values, as 



e.g. those with an Hamiltonian commuting with T[j In this case, in Eq. (10 1 the 
product of operators Y[ G\ X does not depend on the string if we have an even number of 
Majorana operators (i.e. an even number of a x,v ) in each block and it does depend if 
in each block we have an odd number of Majorana fermions (odd-even configurations 
give an odd total number of Majorana fermions that have vanishing expectation value 
by hypothesis). Thus, the density matrix splits in two sums extended only over even 
and odd number of Majorana fermions in each block, i.e. 



1 



PA 1 uA 2 



2 ti- 



[E<n<n<> n<n 

even 



t 



m=o,i leAi ieA 2 \ieAt ieA 2 ! 

odd 

2^ [£{0i0 2 )0l0\ + Y,(OiS0 2 )OtSO\] , (21) 

" even odd 



where S is the complete string of a z belonging to the region between the blocks (19). 
We repeat that the two sums are intended over all possible products of Majorana 
fermions belonging to each interval and even/odd refers to the number of /ii = 1, i.e. 
signaling the presence of a Majorana operator. We introduced also the 'short' 0\ t i for 



a general product of Majorana operators belonging to A\^- Eq. (20) is an immediate 
consequence of the fermionic representation of spin variables and of the considerations 
above: the string S appears only in operators with odd number of each block. 

Spin variables on different sites commute, and so [S, Oi^] — 0. Thus the spin 
RDM can be written as 

l+S l-S 

PA t uA 2 = ~^—P+ H 2~ p - ' ( > 

where p± are the fermionic RDMs 

p± = 7^r 2 [J2(Oi0 2 )olo\ ± J2(OiS0 2 )olo\] = Pcvcn ± Podd . (23) 

even odd 

We will be interested in the Renyi entropies and in writing them in terms of the 
matrices p±. The orthogonal projectors (1 ± S)/2 commute with p± and so the n-th 
power of the spin RDM is the following combination of the n-th powers of fermionic 
RDMs 

Pmua 2 = ^Pl + l -^P- ■ (24) 
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We note that p+ is unitary equivalent to p_, (the complete string Si (or S2) of 
o~ z belonging to Ai (A 2 ) is the unitary operator mapping one into the other, i.e. 
Sip+Si = p-). As a consequence we have 

TrpXuA 2 =^Pl, (25) 

where we used Tr S = 0. Still the matrices p± are not of the form ( 14 1 that are the 
only objects we are able to deal with. Some further manipulations are still needed to 
bring them in a useful form. 

However, before continuing on the main road, we notice that p± are the RDMs 
in fermionic variables of the state 

|*±> = (^ + ^S 2 W >, (26) 



2 2 

where |^o) is the ground state. These are the ground states of the Hamiltonians 

fi±s its \ TT a±s , its 

U ± = n 1 n ^2 



We remark this unitary transformation is nontrivial only for operators on the contact 
surface between the blocks and the remaining chain. Summarizing we can write a 
fermionic RDM equivalent to the spin RDM at the price of adding a finite number of 
non local terms to the Hamiltonian H. 

Now we write the fermionic RDM p± in Eq. ( 23 1 as linear combinations of 
quadratic RDMs of the form (231 that are the ones we are able to deal with. The 
basic intuitive reason why this is indeed possible is that the a z string is the exponential 
of a quadratic form 

Y[<rf =n( 2c i Q_1 ) =exp(wrX) c * c =exp (^^- 1<l2j ' ^ 
11 1 j 

and, in the case of a quadratic Hamiltonian, it can be interpreted as a term added to 
the exponential factor W nm in Eq. (14 1. These observations show that Wick theorem 
applies in some form for Eq. (231, but the 2-point correlations of p dd are different 
from the 2-point correlations of p e ven- 

The RDM's part p dd includes multipoint correlations in which the string S is 
inserted. To deal properly with this kind of objects we introduce the operator (that 
can be thought as a fake density matrix) 

4 = ™ffl, (29) 

where | V I / )( V 1'| is the density matrix of the full system and the denominator is added to 
ensure the normalization Tr^ p^ = 1. This fake density matrix could be introduced 
for any operator, but in the following we need only to consider the string S. If O is an 
operator acting on A, p^ satisfies the property Tya[p^O] = (OS)/(S). If the operator 
S commutes with any operator with support in A, then p^ is hermitian. 

In Eq. (23) the sums are done only on even and odd numbers of Majorana 
operators in each block. We need to write these sums in terms of all Majorana 
fcrmions. A useful observation is that if we change sign to all the a; with I € Ai 
(or I £ A2) in the multipoint correlators, then Oi (O2) changes sign in the odd sum, 
but not in the even one. By taking the appropriate combination of them we are left 
with a sum over all possible Majorana fermions and not only over the even/odd ones. 
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The operator that makes this useful change of sign is the Jordan- Wigner string Si 
(S2) restricted to the first (second) interval, in fact using also S^ 1 = Si, we have 

{—ai , I £ Ai , , v 

, . , 30 

Thus we finally arrive to (1 is the identity matrix) 

_ Pa + S 1 p\Si _ p s A — Sip s A Si 

Pcvon — 2 ' Podd — ^ ' (."*■) 

and so 

Pa + s iPa s i , Pa ~ SipjSi ,_„v 

P± = 2 ± 2 • (32) 

Notice that p\ and Sip\S\ have the same spectrum but different eigenvectors related 
by the matrix S\. The same is true for the matrices p\ and S\p s A S\, but they can 
have eigenvalues smaller than and so in no way can be seen as true density matrices. 
In I Appendix B|we generalized this form to the case of more disjoint intervals. 



Eq. (32 1 is the main result of this section: the rewriting of spin RDM as a 



linear combination of four fermionic RDM, i.e. exponential of a quadratic form as in 



Eq. (14 1. These four matrices do not commute and so they cannot be diagonalized 
simultaneously to find all eigenvalues of the spin RDM. However, if we are interested 
in Renyi entropies with integer a, we can handle this problem in a constructive way: 
we determine the product rules between RDMs and then we construct recursively any 
finite order Renyi entropy. Note that the RDM of the free-fermions is just given by 



p\ in Eq. ( 32 1 . This matrix can be simply diagonalized in the same way explained 



above for the single interval as already done in Ref. [8]. 



2.1. Product rule 



In this subsection we analyze the algebra of RDMs generated by a quadratic form, i.e 

1 

Z{W) 



Pw = yfw\ exp (5Z a i w ^ a n 1 4 ) > (33) 



Ln 



where we do not assume W to be hermitian (to include the non-hermitian contribution 
of the string). Anti-commutation relations of Majorana operators make always W a 
complex skew-symmetric matrix, i.e. W T — —W. The constant 

Z(W) = Trexp(j] ai W ln a n / 4) , (34) 

Ln 



ensures the normalization Trp^ = 1. In Appendix C we show that this normalization 



for complex diagonalizable skew-symmetric matrices is 

Z(W)= [] 2cosh(^), (35) 



where {w}/± is the set of eigenvalues of W with halved degeneration (W is a skew- 
symmetric matrix so any even function of W has eigenvalues with even degeneracy). 



We assumed Z(W) 7^ 0, pathological cases can be cured as explained in Appendix D 



The product of fermionic RDMs of the form ( 33 1 is 



Z(\og(exp(W)exp(W'))) 
Pwpw - Z(W)Z(W) Piog(cxp(w)cx P (W')) ■ (3o) 
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Indeed the commutator of operators in the exponent of ( 33 ) 
V- W ^ W 'ik L . 1 _ a T {W,W'}a 



^ 16 

l,n,j,k 



aia ni djdk 



(37) 



is the essential ingredient in the Baker-Campbell-Hausdorff formula for the product 
of exponential of operators. 

Fermionic RDMs are specified by correlation matrices 

Tij = Tr[(npwaj] - , (38) 



which can be written as (see Appendix C ) 

r = tanh(f) =* e- = ^. (39) 

Clearly the second equation is true only when 1 — T is an invertible matrix. At this 
point, let us briefly summarize the logic of the following derivation. We can easily 
calculate/manipulate the correlation matrix T that via Eq. (39) gives the exponential 
factor W that defines the quadratic density matrix in Eq. ( 33 ) . We need to find what 
are the consequences of the product rule of RDMs for the correlation matrices, i.e. 
we need to find what is the correlation matrix corresponding to the product of two 
RDMs. While, through the chains of equations above, any W defines a single pa, the 
opposite is not true and there are several possible W's for each pa- Nevertheless, we 
can give a unique recipe for the composition of correlation matrices. 

We indicate this matrix operation as T x T' (notice it is not the product of the 
matrices) and it is formally defined by Eq. ( 36 1 as 

p[r]p[r'} = Tr[p[r}p[r'}]p[rxr'}. (40) 

To specify this operation we still need two ingredients: 

(i) an usable expression for the correlation matrix 

associated to the product pwPw = /°[r]/?[r']; 

(ii) an expression for the trace of two fermionic RDMs 

fr r'x - t rn rr'i Z(W)Z(W) 

{E, r } = Tr p[T]p[T ] = z{log(eMw)eMwi))) , (42) 

in terms of the correlation matrices T and T' . 

The first requirement is easily obtained if we assume 1 — T and 1 — V invertible. 
Indeed, if we make explicit the exponential products 

i + rxr _ i + n + r 

i-rxF _ i-n-r'' ^ ' 

after simple algebra we obtain 

rxr' = i-(i-r') T ^ ff7 (i-r). (44) 

r x r' is a skew-symmetric matrix, even if it is not obvious from the above formula. 
We checked that this relation remains true also if 1 — T is not invertible (at least for 
the kind of matrices we are interested in) , but a complete rigorous proof of Eq. ( 44 1 
is beyond the goal of this paper. 
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The second request is less trivial because the correlation matrix T does not 
determine univocally the matrix W, and the sign of Z{W) remains ambiguous. 
However, {r, T'} is a functional of T and V , i.e. it is the product of the eigenvalues of 
(1 + rr')/2 with halved degeneration (the spectrum of IT" is double degenerate [2*5]) 



{r,r'} = 



n 



i+fi 



^eSpcctrum[IT']/2 



= ±A/det 



1 + TF 



(45) 



The unspecified ± sign in front is the ambiguity that (as shown in Appendix C ) can 
be solved by rewriting the composition rule as 



{r,r'} = ex P (i 



dA IT' - 1 
Tr 



1 + A AIT' + 1 



(46) 



'7[o-+i] 

that does not depend on the curve 7. 

It is evident that the operation x is associative and so we are in position to make 
any product of fermionic RDMs: 



jjp[r i ] = {r 1 ,...,r n }p[r 1 



x r„] 



(47) 



where 

{r x , r 2 , r 3 . . . , r n } = Tr^Mr,] ■■■} = {r x , r 2 }{r x x r 2 , r 3 . . . , r„} . 



(48) 



If, for some i and j, the matrix (1 + FiTj)/2 is not invertible, Eq. (47 1 cannot be 
applied. It is still possible to isolate the pathological parts and we report the details 
of this procedure in |Appendix D[ 

3. Renyi entropies 



From equations (25), (32) and (47), the Renyi entropies for integer a can be written 

as 



s = — l0 4i E n c [C,]{r Cl ,---,rcJ 



1 - a 



(49) 



where we defined the variables Q = 1> 2, 3, 4 (that label which of the terms in Eq. (49 ) 
is taken in the particular product) and we defined the shorts for the 2-point correlation 
matrices 



ri = i> 

and 



r 2 — r Slp i Sl , 
4Q = { (s) 

-(S) 



r 3 = r p s , 



Ce {1,2}, 

C = 3, 
C = 4. 



r 4 = r 



Sip s Si 



(50) 



(51) 



In the case of a — 2 there are no problems with singular terms (see Appendix C ) 
and the above expression can be re-written as 



S 2 = - log 



^cICrMCsKiWdet 



Cl,<2 



1 - r Cl 1 - r 



C2 



1 + r Cl 1 + r C2 



(52) 
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where here and in the following equation we leave the sign ambiguity unspecified. 
Taking into account the trace's invariance under cyclic permutations, S2 becomes the 
logarithm of a sum of 10 terms 



<dct 



1 



\ £ C [CiMC 2 ](±)Jdet 



Cl>C2 



1 - 17, 1 - r 



C2 



1 + r Cl 1 + r <2 



(53) 



Pi 



(54) 



These formulas are already usable for a direct computation of Renyi entropies. 
There are however some simplifications that occur by using the property of the 
correlation matrices Vs. The first matrix 1^ is the standard fermionic correlation 
matrix (i.e. the one corresponding to free fermions in the absence of the Jordan- Wigner 
string, already considered in Ref. [5]), r 2 can be obtained from r 4 as r 2 = P\TiPi, 
with P\ the Hermitian unitary matrix 

u 2 

The same relation occurs between the third and the forth matrix r 4 = P 1 r 3 P 1 . Instead 
r 3 is not trivially related to T\. In |Appendix A| we prove the following identity 

r 3 = Tx - T ABl v B 1 iB r BlA , (55) 

where the double subscripts take into account restrictions to rectangular correlation 
matrices, i.e. the first (second) subscript identifies the region where the row (column) 
index runs. We show in |Appcudix B| that similar properties are valid also for an 
arbitrary number of intervals. 

Using these relations after some algebraic manipulations one can write down the 
full sums for the Renyi entropies. It is important to notice that (S) = Pf(rs 1 ) and so 
(S) 2 = det(T Bl ). Furthermore to short the notations we denote 

{ • • • • I !'• • • • 1 {•••.! ,•••.!,.•••! . (56) 



Finally S2 can be written in a rather simple way 

■{r?} + {r l5 r 2 } 



s 2 



log 



+ (det r Bl ) 



{r|}-{r 3 ,r 4 } 



(57) 



2 2 
But, increasing the order a, the explicit expressions become soon long. For example, 
we give the (simplified) formulae for S3 



S3 = 
and Sa 



S 4 = - - log 



+ 3(detr Sl )- 



{r 1 ,rl} + {r 2l rl}-2{r 1 ,r 4 ,r 3 } 



{r3} + 3 {r?,r 2 } 

4 ■ -v— -1, 4 

{rf } + 4{r?, r 2 } + 2{r?, r^} + {r x , r 2 , r x , r 2 } 



(detr Bl ) 



8 

/ {r 1 ,r 3 ,r 1 ,r 3 } + {r 1 ,r4,r 1 ,r 4 } {r^,r§} + {rj,rl} 
I 4 2 



{r^ra^.^g} - {r 3 ,r 1 ,r 4 ,r 1 } - {r 1 ,r 2) r 3 ,r 4 } 



(58) 
(59) 



-{r 1) r3,r 2 ,r4}-{r 1 ,r 2 ,r4,r 3 } 



{r 1 ,r 2 ,r2}-{r 3 ,r 4 ,r2} 



(detrs )2 {r|} - 4{r|, r 4 } + 2{r|, r|} + {r 3 , r 4 , r 3 , r 4 } 
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4. Critical XX model 

In this section we report the explicit results for the XX chain in zero magnetic field (i.e. 
7 = and h = in Eq. ([9])). XX chains have been previously analyzed by Furukawa et 
al [TT] by means of exact diagonalization techniques that allowed to explore relatively 
small chains with at most 30 spins. The asymptotic results from CFT (cf. Eq. ([5])) 
were obscured in this previous analysis by large oscillating corrections to the scaling. 
The smallness of the systems and the lack of a precise knowledge of the form of the 
correction to the scaling made impossible any finite-size scaling analysis to check the 
CFT predictions In fact these oscillating corrections to the scaling have been 
widely observed in last years [313 123 I2E| , but a precise theory about their origin and 
their exact form is only recently available [2"4"1 |2"9"] . 

By exploiting the exact solution, we can avoid these problems and explore large 
enough values of £ allowing a finite-^ scaling analysis similar to the one for a single 
block [24] . We start from the infinite volume limit. We consider two intervals both 
of length £ — l\ = £2 at distance r. In this case, the four-point ratio x in Eq. Q is 
given by 

2 

(60) 



J + 

and the scaling of the quantity determining the Renyi entropies is (we use c = 1) 

/ (£ + r\ 2 \ fa" 1 /")/ 6 

Tr ^ c H^T7)) (61) 

that is an implicit definition of F^ t (x,£) encoding the fact that for finite £ we expect 
corrections to the scaling of the form 

F^(x,£) = F£ FT (x) + £- s °f a (x) + . . . , (62) 

where F^ FT (x) is the universal quantity appearing in Eq. The exponent S a 

governs the leading correction to the scaling. It has been shown with CFT in Ref. [29] 
that this exponent is equal to 2x/a independently of the number of intervals and x is 
the scaling dimension of an operator introduced by the conical singularity necessary to 
describe the reduced density matrix. For a single interval, it has been shown by exact 
analytic calculation that S a — 2/ a |24j and thus we expect the same exponent for the 
double interval case (the presence of a relevant operator with dimension x — 1 has 
been justified in Ref. [29]). In the following analysis we will greatly benefit of the exact 
knowledge of the quantity c a from the exact solution of the single interval entanglement 
[2"51 15U] and so our results will be also a non-trivial check of the equivalence of the two 
non-universal multiplicative constants in the single and double interval case. 

We report in Fig. [^results for the function Fl a ^ 4 (x,£) for various £ and x. It is 
evident that irrespective of the value of x, with increasing £ the results approach the 
CFT prediction. For odd £ the asymptotic result is approached from below, while for 
even ones it is approached from above. These are the already mentioned oscillations 
that made difficult the analysis based on small chains. These plots are the first direct 
quantitative tests of the correctness of Eq. (jfjj). These figures do not leave doubts 
about the correctness of the results of Ref. [12]. It is worth to mention that, as 
already shown elsewhere [Til E2J H3] , the finite £ curves do not have the symmetry 
x — y 1 — x that is restored only in the £ — > 00 limit Eq. Q . 

Having established the correctness of the asymptotic form, we can now move to 
the finite £ corrections and check whether the prediction S a = 2/a [55] is correct. 
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Figure 2. Scaling function F^ at (x) for a = 2, 3, 4 (from top to bottom) for the XX 
model in the thermodynamic limit and for various £ = i\ = I2 and r. Corrections 
to the scaling show even-odd oscillations with I as in the single interval case. The 
results converge quickly to the universal CFT prediction F a (x). 
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Figure 3. Corrections to the scaling for F^ l (x) obtained by subtracting the 
asymptotic value. The first three panels show that for a = 2, 3, 4 at fixed x = 1/9, 
both for fermionic and spin variables corrections to the scaling have exponents 
6 a = 2/ a, as for one interval. For fermionic degrees of freedom the asymptotic 
behavior is obtained for smaller values of I. Last panel: x independence of the 
corrections to the scaling exponents. All curves with different x at fixed a lie on 
parallel lines. 

Having precise control of the corrections to the scaling is not an academic task: their 
analysis is fundamental to provide accurate results when such large system sizes are 
not available (in previous studies they would have been an important tool) and in 
cases when the asymptotic form is not known (as in the following section for the 
Ising model for a > 2). In Fig. [3] we report the function \F a {x) — F^ t (x)\ for fixed 
x = 1/9 and a = 2,3,4. We both report results for the spin Renyi entropies and 
for the fermionic ones (i.e. without considering the string contribution). We recall 
that for free fermions we have F a (x) = 1 identically. The results show a power law 
behavior for large enough £ with the predicted exponent 2/ a as for the single interval. 
Notice that by increasing a, the values £ where the leading asymptotic correction can 
be identified become larger and larger, in analogy with the single block case [53] (as 
obvious because of the smallness of the exponent S a ). For the fermionic variables the 
asymptotic behavior is reached before than for spin degrees of freedom: the string 
introduces further corrections to the scaling that in the present model are subleading. 
To show the x independence of this exponent in the last panel of Fig. [3j we report the 
same kind of plots for different values of x, showing that, at fixed a, the corrections 
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Figure 4. The function Fj (a;) for a small chain with L = 39 spins. Oscillating 
corrections to the scaling prevent to see the universal CFT prediction shown as a 
continuous curve. 



lie on parallel lines. 

To conclude this section we present some results for finite systems. It is known 
that all the formulas above (including the x dependence) at the leading order in 
finite systems can be described by replacing any distance Uij by the chord length 
Uij — > ^ sin We checked that this rescaling is indeed correct and that, for large 
enough £ and L, the results agree with CFT. These plots give no further information 
compared to the ones already presented and we do not report them. We only show 
the results for a rather small system of length L = 39 (that nevertheless is above 
anything obtainable by exact diagonalization) . In Fig [4] we report the resulting 
for all the possible divisions in four parts of this chain of length L = 39. For such 
small chain, the results are obviously very unclear since the corrections to the scaling 
are obscuring the CFT scaling represented by a continuous line that is surrounded 
by the points, signaling the oscillatory nature of the corrections (for clarity, compare 
with the analogous plot for the Ising model in next section). 

5. Critical Ising universality class 

The Ising model is given by Hamiltonian ([9]) with 7=1 and it is critical for h = 1. 
Results for this model have been already derived numerically for a = 2 in Ref. 
[13] by using a tree tensor network algorithm and Monte Carlo simulations of the 
two-dimensional classical problem in the same universality class. These calculations 
allowed a precise determination of F2(x), but the system sizes explored were not 
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Figure 5. Scaling function F^ l (x) for a = 2,3,4 (from top to bottom) for 
the Ising model in the thermodynamic limit and for various £. Corrections to 
the scaling are monotonous. The top curve in each plot is the extrapolation to 
I — > oo. The convergence to the universal CFT prediction F a (x) is slower than 
in the XX case, because the leading exponent of corrections to the scaling is 1/a. 




Figure 6. F^ orm (a:) for different x and 7. All results present the same leading 
correction to the scaling exponent 82 = 1. The extrapolated data at I — > 00 
collapse in the single point F2(x) = 1. 



enough to analyze Renyi entropies with larger values of a. In the course of our 
analysis we always compared our data for small systems with those in Ref. |13j . in 
order to check the correctness of both methods. 

We report our results in the thermodynamic limit for various values of I = i\ = £2 
at different separations r (resulting in the four-point ratio x given in Eq. (60)) in Fig. 
[5] As already noticed in Ref. [T3|, oppositely to the XX chain, in the Ising model we 
have monotonic finite I correction to the scaling. For finite £, the results do not show 
the symmetry x — > 1 — x valid for infinite I. This is restored only by the extrapolated 
data at I —¥ 00. To perform this extrapolation in the most accurate way, we have first 
to determine the correction to the scaling exponent. For a single interval, it is exactly 
known that the leading corrections are characterized by the exponent S a = 2/ a as for 
the XX chain (the origin of the equality [24] is due to the equivalence of the two RDMs 
[31]). However for a = 2, it has been already shown [13] that 82 = 1/2, different from 
the single interval one. This result is quite surprising also in view the recent CFT 
analysis |29| predicting the same behavior for any number of intervals. Thus as a 
first step we check the corrections to the scaling exponent for the Ising model in the 
absence of the Jordan- Wigner string between the two blocks (i.e. we consider only 
the correlation matrix Ti, as done in Ref. [8]). In this case, the results reported 
in Fig. [6] give a compelling evidence that the leading corrections to the scaling are 
given by S a = 2/ a as in the single interval (data not reported here for a = 3,4 show 
that this remains true). At this point it is natural that the Jordan- Wigner string 
produces another operator at the conical singularity, that in the Ising model is the 
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Figure 7. Universality with respect to 7. We plot Fl at (x) for x = 1/4 and 
x = 3/4 for different values of 7 in the Hamiltonian. All results present the same 
leading correction to the scaling exponent 82 = 1/2. The extrapolated data at 
£ — > 00 collapse in a single point equal to F 2 (l/4) = F 2 (3/4) = 1.31886 . . .. 



leading one. According to Ref. all the corrections to the scaling should be of 
the form 6 a = 2x/a, thus taking the result 62 = 1/2 for granted (see also below for 
further evidences) , we conclude that the Jordan- Wigner string introduces an operator 
with scaling dimension x = 1/2. Such an operator in the continuum limit of the Ising 
model exists and it is the Majorana fermion, that has exactly the same features of 
the Jordan- Wigner string (i.e. same symmetry and same non-local character). Such 
an operator is clearly not present in the single interval case. These considerations 
allow us to conclude that the leading corrections to the scaling for the double interval 
entanglement in the Ising model are described by the exponent 

5 a = -. (63) 
a 

Unfortunately already for a = 3, the value of £3 = 1/3 is very low and subleading 
corrections to the scaling going with exponents mS a (with m integer) are expected 
to influence the results in a considerable manner. For this reason, in order to have 
an accurate determination of the asymptotic behavior, at fixed x we consider all 
corrections to the scaling up to those with exponent 1. The resulting extrapolated 
data are the top points in Fig. [5J In the case of a = 2, the CFT prediction F%{x) (cf. 
Eq. perfectly with the extrapolated data agree, giving strong support both for 

the procedure to account for the subleading corrections terms and on the asymptotic 
form. For a 7^ 2 no CFT prediction is still available. The extrapolated data in Fig. [5] 
are the first data for infinite t. 
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We now turn to consider the issue of universality. All the critical models (h = 1) 
for any value of 7 7^ are in the Ising universality class. However, the results at finite £ 
show a strong dependency on 7 (as obvious from the different correlation matrices) . In 
Fig. [7] we report several data for F^^x) for different values of 7 and I at fixed x = 1/4 
and x = 3/4 (that for the symmetry x — > 1 — x have the same asymptotic value). At 
finite £, all results are evidently different. In the figure we report the extrapolation 
with two corrections to the scaling (i.e. with 62 = 1/2 and 262)- For £ — > 00 all data 
tend to the same value predicted by CFT F 2 (l/4) = (2 + \/2(l + 3 1/4 ))/4, confirming 
in a single plot many results: (i) universality with respect to 7, (ii) correctness of the 
correction to the scaling form, (iii) correctness of the CFT prediction Eq. Q . 

To conclude this section we report the data for a finite chain. As usual, in all 
the scaling variables we substitute distances with the chord distances. In Fig. [HJ we 
report the values of F2 at (a;) for all the possible choices of intervals A± t 2 and B12 in a 
chain of L = 39 spins. Compared to the analogous plot for the XX chain (Fig. [4]), the 
figure is much clearer, due to the fact that corrections to the scaling are monotonous. 
Notice however that the data points lie much below the asymptotic value because the 
exponent 82 = 1/2 is small (also compared to the XX case with 82 = 1) resulting in 
very large corrections to the scaling. 




Figure 9. Left: Comparison of the Renyi entropy 52 for a non-critical Ising 
model of two disjoint intervals versus the double of the single interval and the 
fermionic one. While the fermionic entropy is asymptotically always the double 
of the single interval case, the one in spin variables presents a — log 2 difference due 
to the string. Right: Mutual entropy AS2 as a function of the magnetic field h for 
an Ising chain of 129 and 257 spins. The crossing at the phase transition is shown 
for three configurations. The intercept with h = 1 is in good agreement with 
the CFT prediction. The opaque lines are the corresponding fermionic mutual 
entropies 



6. Non critical models 

The richness of the phase diagram of the XY model allows us also to explore gapped 
phases that are almost everywhere except on the line \h\ = 1 and the segment 7 = 
with ft 2 < 1. These systems have not be considered so much in the literature until 
now, because for non-critical models all the correlations and entanglement between 
the two blocks fall off exponentially (with a decay rate given by the inverse gap or 
mass). Thus one would always expect 

S a (£, r,l') = S a (£) + Sail') + 0{e~ r ' A ) . (64) 

However, this is not so obvious because of the importance of the connected part in 
the correlations. In Fig. [9] (left) we report S2 for the double interval case. While for 
ft, > 1, the spin entropy is the same as the double of the single interval and the same as 
the fermionic one, for ft, < 1 there is clearly an offset (that we quantify in — log 2, see 
below) . The best way of detecting these unexpected effects is to consider the mutual 
entropy 

A$ Q (4r,0 = S a {i) + S a (f) - S a (£,rJ') , (65) 

that gives automatically zero when factorization occurs. 

From numerical data we deduce that, for large blocks, all terms 

a 

nc[Ci]{r ClJ ...,r Co }, (66) 



in Eq. (49) have the same absolute value. However, they can have different signs. 



Furthermore numerical evidence suggests that terms consisting only of the correlation 
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matrices Ti and T2 have always positive sign. II these observations are generally true, 
we have 



AS a l — log 



1 + ^ E ("i^KC}]], (67) 

{C} S.t. 
#3+#4^0 

where e[{C}] is the sign associated to the element {T^, • • • , T^} and #3(4) is the 
number of correlation matrices r 3 ( 4 ). We expect an eventual discontinuity in AS a 
when crossing a critical line. Thus we study non-critical chains with magnetic field 
close to h = 1. We found numerically that only terms with odd number of correlation 
matrices ^(4) display sign changes. Thus (see right panel of Fig. [9] for the explicit 
plot) we conclude from the numerical evidence, the behavior 

flog2 h?<l, 
AS a = \ B (68) 
I otherwise . 

This result is not completely unexpected: also Sa of the single interval for \h\ = 
tends to log 2, independently on £ [5] [231 [3D] . Thus in the definition of AS a , since also 
the double interval S a tends to the same value, we are left with a single log 2j|] 

We stress that at h = 1 the various A Q cross in a single point if they are 
characterized by the same four point ratio x (see right panel of Fig. [9]) that is the CFT 
prediction. These kind of plots could be used to detect the phase transition points in 
systems where it is not exactly known. 

The a independence of the previous expression, allows us to analytical continue 
the result to a = 1 and to conjecture the same behavior for the asymptotic von 
Neumann entanglement entropy in non-critical regions 

Si = Sf Im - 0(1 - h 2 ) In 2 + . . . , (69) 

where S^ Tm is the fermionic entanglement entropy obtained from the correlation 
matrix T\. Notice that this does not exclude that exponential corrections to this 
asymptotic form in the fermionic and spin variables could have different amplitudes. 
It would be interesting to explore this issue with the continuum theory in the form- 
factor approach [52"] . 



7. Entanglement evolution following a quench 

Finally we consider the evolution of Renyi entropies after global quenches. In this 
problem, the system is prepared in the ground-state l^o) of an Hamiltonian Hq of 
the form ^ and then it is let evolve from time t — with a different Hamiltonian 
H, always of the form ^ but with (ft., 7) ^ (/io,7o)- The main feature of this non- 
equilibrium problem is that the initial state differs globally from the ground state and 
the excess of energy (compared to the ground-state of H) is extensive. A connected 
spin block reacts to the quench increasing Renyi entropies linearly in time up to 
each spin in the subsystem becomes entangled with the environment. Then entropies 
saturate. There is no reason why disjoint blocks should behave in different way, except 

X This is a consequence of the double degeneration of the ground state for \h\ < 1. It is easily 
understood at h = 0, where the ground state is any linear combination of the states all up and 
all down that we can denote with | j) and | -l). Since there is no symmetry breaking term in the 
Hamiltonian, the diagonalization selects a state with zero magnetization, i.e. (| t) ± | 4,})/v2, in 
which the entanglement of any subsystem, connected or not, is always log 2 as stated in the main 
text. 
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for the necessity to subtract the mutual entropy. Our only goal here is to understand 
the differences between the spin representation and the fermionic one, postponing 
any more accurate analysis to further studies. We only report results for the Renyi 
entropy S2, but the value of a is unimportant. We compare the numerical data with 
the prediction that follows from the interpretation of the entanglement evolution in 
terms of motion of quasiparticles [33] . According to this physical scenario, because the 
excess of energy of the system is extensive, any site acts as a source of quasiparticle 
excitations. Particles emitted from different points (further apart than the correlation 
length in the initial state) are incoherent, but pairs of particles moving to the left or 
right from a given point are entangled. Thus S a (t) should just be proportional to the 
number of coherent particles that emitted from any point reach one a point in A and 
the other a point in B. Since there is a maximum speed for these excitations %, for 
a single interval this implies the linear growth for 2vjMt < £ and saturation for very 
large times. For a general bipartition the result is [33J 

S a (t)m [ dx' { dx" { dx f H a (p)dp6(x'-x-v(p)t)S(x"-x-v(p)t) ,(70) 

Jx'eA Jx"£B J-oo J 

where we assumed momentum conservation. For a double interval, this formula 
predicts a series of linear behavior with different slopes, that finally saturates at late 
time [33]. The function H a does not depend on the subsystem length and topology. 
It has been exactly derived for a single interval [34] 



H a (cos A) 



1 



1 — Q 



log 



1 + cos A \ a 



1 — cos An a 



(71) 



2 / V 2 

where A(p) is the difference between the Bogolioubov angles before and after the 
quench, whose explicit expression in term of the dispersion relation e(ip) is reported 
in |34j . Thus, plugging everything together, the leading order of Renyi entropies is 

S iur ,e 2 {t) ~ s tl (t) + S t3 (t) - AS[h,r,mt), (72) 

with (x[ a ,b] ( x ) = d{x - a)0(b - a;)) 



AS(t) 



and 



v dip r 2 K\* , 

— fla(cOsA) / dv^X[r,r+min(£i,t2)]( v ) ~ X[r+max(ii ,t 2 ) ,r+ti+£ 2 ] ( v ) 



S e (t)=t[ ^2\e'\H a (cosA) +1 [ JtJ q (cosA). (73) 

J2\e'\t<t 2?r J2\e'\t>e 27r 

Let us now go back to the direct computation of the entanglement in two blocks 
after a quench, we observe that not all terms contribute to the leading order in blocks 
lengths and distance: the string expectation value decays exponentially and only terms 
constructed with T\ and T2 survive. For instance the Renyi entropy S2 is simply the 
logarithm of two terms 

r x iY 



5'- 



log 



1 /l , 1 
-exp(-Trlog- 



>logi 2 



1 f\ , 1 

-exp^-Trlog — 



r? 



1 Trlog- 

2 B 2 



rir 2 



log 2 + c 



+ c + log 2 



if Trlog iiTiIa < Trlog Id^l 

if Trlog 1+r 2 lT2 - Trlog 
otherwise. 



(74) 



(75) 
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Figure 10. Evolution of the Renyi entropy 5*2 after a global quench from a 
non critic al s ystem to a critical Ising. The asymptotic prediction is given by 
equation ( |72| while the fermionic entanglement is obtained neglecting the string 
contribution. When {Ti,^} becomes negligible with respect to {T^}, a jump 
of about ~ log 2 is seen. In the graph this happen in the neighborhood of the 
time t ~ 11. The asymptotic prediction presents a small offset due to the finite 
entanglement in the initial state. 



Analyzing several data (we report only a single example in Fig. 10 1, we deduce that in 
the scaling limit the term constructed with the standard fermionic correlations is never 
negligible: finite order Renyi entropies are controlled by fermionic correlations. When 
other terms are comparable with {Ti, ■ • • ,Ti}, they give additive O(l) contribution 
of the form ~ logfc, where k is the coefficient in front of the factor {• • •}, and k = 2 
for S 2 - 

The quasi-particle interpretation |33j perfectly agrees (up to a constant) with the 



fermionic representation (see Fig. 10). In other words the non locality of the 



Jordan- Wigner transformation does not influence (as one expects) the leading order 
of time evolution of entropies after a global quench. However, as in the non-critical 
equilibrium case, we observe an additive log 2 contribution (in Fig. 10 a stick of width 
log 2 is included to appreciate this difference) from the string and we believe that this 
is a general feature that will persist also away from criticality. As a final comment, 
we mention that the O(l) offset of the numerical data at finite £ compared with the 
asymptotic form is a consequence of the initial entanglement. This subleading term 
is present also for the single interval case [31], and can be included in the asymptotic 
result 135 . 
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8. Discussions 



We presented a general method to calculate Renyi entropies for integer a in the ground 
state of the XY Hamiltonian ^ in the case when A consists of two (but the method 
works for even more) disjoint intervals. We carefully analyze the results in the critical 
cases of XX and Ising universality classes. We found that the asymptotic results 
for large intervals and separations are described correctly by conformal field theory 
when results are known. We also provided results for the asymptotic scaling of the 
a = 3, 4 Renyi entropies of the Ising model, for which no analytic prediction is still 
available. In order to check CFT predictions we had to carefully include corrections to 
the scaling in the analysis. We found that for the XX model, the leading corrections to 
the scaling are governed by the exponent S a — 2/ a, the same as the single interval case 
[24] , in agreement with general predictions [29] . Oppositely for the critical Ising model 
we found S a — that is half of the single interval one. This can be understood 
in terms of the general theory of the corrections to the scaling [53] by interpreting 
the Jordan- Wigner string in the lattice model as a non-local fermion in the Ising 
CFT. We also have briefly considered the entanglement in the gapped phase and the 
time evolution following a quantum quench. The provided exact representation of the 
reduced density matrix of spin variables could also be useful for the calculation of the 
entanglement between Ai and A 2 measured by the negativity [15] . 

A major problem remains still open. We are unable to find all the eigenvalues 
of the reduced density matrix of two intervals and so the Renyi entropies for non- 
integer a and in particular for a — 1 that would give the widely studied von Neumann 
entropy. The same problem is also present in conformal calculations. It is possible 
to obtain Tr p\ only for a integer [T2] and the analytic continuation of the result to 
general complex values remains a big open problem. 

The contribution of the Jordan- Wigner string also affects the entanglement 
entropy of a single interval in a system with one or more boundaries, if the block 
A does not include the boundary. In this case, the correlation matrix is slightly more 
complicated and we are currently analyzing the problem. Among the other things, 
this would be relevant also for some quench problems where indeed only the fermionic 
entanglement has been calculated [31)] . 

Another interesting open question concerns the behavior of the Renyi entropies 
in systems with disorder for which in the single block case it is known that the 
entanglement simply follows from the counting of the singlets [37] between A and 
B. Whether this remains true in the case of more intervals is not yet known. 



Appendix A. Correlation matrices 



We report in this appendix some simple properties of the correlation matrices 

If - = Tt[Oi Ail - i, j € A , (A.l) 



with ps defined in Eq. (29 1. r s 's elements are ratio of two expectation values (n ^ 0) 

ai+n UjeBi ( a 2j-io,2j)ai 



i ll+n — 1 \ i (A-2) 

which can be evaluated by means of the Wick theorem. We show that V s can be 
expressed in terms of the system correlation matrix T = r , where 1 is the identity 
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matrix. By isolating the two-point correlations between a/ and the other Majorana 
operators in the numerator of ( |A.2[ ) we obtain 



a i+n \\ {aaj-ia2j)ai\ = (S)T u+n + [TAT] 



jeBi 



J l l+n ' 



with 



A„ = (-l)^sign(j-0( 1] 



keB' 



(A.3) 



(A.4) 



where B*j denotes the set of all indices {21} U {21 — 1} with I G Si and with i and 
j removed. The expectation value in (A.4) is the Pfaffian of the system correlation 
matrix restricted to the region B*j . Checking the Pfaffian sign we have 



A, 



with Tb ± the correlation matrix restricted to the region B\ 
expression into (A. 2) we get 

T s = T A 



1 Bi 



AuB 



i/ r B! , 



(A.5) 

Substituting this 
(A.6) 



where the double subscript takes into account restrictions to rectangular correlation 
matrices, i.e. the first (second) subscript identifies the region where the row (column) 
index runs. Notice that T s is the Schur complement of in Taub-i ■ This proves 
Eq. pH). 



Appendix B. Generalization to many disjoint spin blocks 

We consider in the appendix how to generalize our approach to the case when the 
subsystem consists of n disjoint blocks of lengths ii, with i = 1, ... ,n. We indicate 
with Of a product of an even number of Majorana operators lying in the i-th block 
Ai and with 0~ an odd product. Sf is the complete a z string associated to the i-th 
region between consecutive blocks Bi and Sf is the string corresponding to the block 
Ai (see Fig. [ljfor the double block case). Eq. (20) can be easily generalized observing 



that the reduced density matrix has the following expansion 



f 



Pu iAi = ^j:Yl (S[{s}]0?---0^)(0?--.O s n «yS[{s}}, (B.l) 

where s, = ±1 and ^fls}] are a z strings determined by the conditions 

fe-l 

S[{s 1: ■ ■ ■ , s k , ■ ■ ■ s n }] = S[{ Sl , • ■ ■ , -s k , ■ ■ ■ s n }] 11 Sf , (B.2) 

i=i 

S[{0,---,0}] =1. 

Indeed, any spin operator that can be written as a product of Majorana operators that 
are odd in number if restricted to the fc-th block is constructed with an odd number of 
fj x (v) f the fc-th block. Spin variables a x and a v are non-local in the fermionic space 

fc-i 

ofJi = A(Ai, . . . , A k ^) l[ Si a 2i(2i _i) , (B.3) 

where A(Ai, . . . , A k _i) is an operator that depends on the Majorana fermions of the 
blocks Ai,... Ak-i- When multiplying an odd number of a x ^ v) the string Jl^i 1 &i does 
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not cancel (we remind Sf — 1) so we obtain the conditions (B.2). The construction 
of the fermionic equivalent reduced density matrix is analogous to the double block 
case. In fact the unitary transformation 



^EIK 



Si G { + , -} . 



with 



n 



i + si 



i-1 



(B.4) 



brings the spin reduced density matrix to a fermionic one. In particular we obtain the 
simple result 



PUiAi 



^"E E (W-O^XOJ'-O?)*- (B.5) 



{a} O'^-O 



The term associated to the configuration {s} 

-^y E (S[{s}}0^---0^)(0^---0^y 

2^ 1 «l... «n 



S[{s}] 

P{s} 



(B, 



is sum of 2™ 1 fermionic RDMs. This is because, fixed 5[{s}], the further dependence 
on Si can be handle in the following way (remind S k O s k S k — sO s k ) 



S[{s}\ 
"{si,---,s fc ,---,s„} 



r SIM] , S[{s}] , 

W'{ai,-,«A„-.*n} + ^{ai>-,-«fc,-,»ft}J 



+ 



lP{ su -, Sk ,-, Sn } + P{ Su ...- Sk ,...,s n }\ A k 



and, for any string S B and S , the operator 



, (B.7) 



(B- 



is a fermionic RDM. We obtain 2™ _1 and not 2 n fermionic RDMs only because <5[{s}] 
does not depend on sj (there is no string contribution for operators lying in the first 
block). Thus we can write 
1 



PUiAi 



E 



pf [ r U 3 , w BjE e «^ r &] **e{-M>. (B.9) 

{z}w/o repetition {s} 

The explicit form of the correlation matrices r|°^ is a simple generalization of the 
double block case 

As} 



r {i} = P I S 1> ■ • • ' S «-l] r .4u(l J B Z )I T \ J B, P [ S 1> • • • > S «~l] 



(B.IO) 



where we indicate with A/B the Schur complement of B in A (see Appendix A). Pf 
denotes the Pfaffian, and P is the diagonal matrix 



P[si, . . . ,S„_i] 



Sn-llA„_i 



(B.ll) 



u. J 
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The index i G {0, n — 1} runs over the environment blocks lying between the first and 
the last subsystem block, and B = so that ^{o} is the standard fermionic correlation 
matrix Ta = Ty A . The sign e is determined by the conditions 

{s} _ f i {0 = {o}v{ s } = {i,...,i}, 

6{l} ~ | | L)* i' 1 ,;;; • ^ otherwise , ' L> ' ' ' 

where (— l)-? e W is equal to —1 if j € {i} and 1 otherwise, and 

{s}sj^— sj = {si, ■ ■ ■ , Sj-i, —Sj,Sj+x, • • • , s„_i} . (B.13) 

Renyi entropies can be computed in the same way as when the subsystem consists of 
just two disjoint blocks. But many more terms contribute. And when the number 
of blocks is comparable with the chain size we expect an extensive behavior of Renyi 
entropies, as observed in [18] , that is reminiscence of the huge number of fermionic 
RDMs needed to represent the spin RDM. Finally, as in the double block case, Renyi 
entropies do not depend on the sign of Pfaffians because any sign change is the result 
of a unitary transformation (that is just a o~ z string). 



Appendix C. Some technical details for the product rules 

The reduced density matrix of a given number of disjoint spin blocks in a free chain is 
unitarily equivalent to a finite sum of operators that well behave under multiplication. 
In fact one of our results is that, for some a z strings Sj and constants Cj, the relation 

Tr A SjPo _ ^ s 



c.^m^c^, (c.i) 



3 J 3 



holds. Like in the single block case each operator p s i turns out to be a normalized 
exponential of a quadratic form in the Majorana fermions. Up to a constant factor, 
these matrices are closed under multiplication and we identify any of them by the 



skew-symmetric matrix W in the exponent in Eq. ( 33 1 . We recall in this appendix 
some basic properties of skew-symmetric diagonalizable matrices. 
If A is an eigenvalue of W, then also —A is because 

= det|W - A/| = det\W T - A/| = det|-W - Xl\ . (C.2) 

If w\ and are two eigenvectors with eigenvalues A and /i, then 

p,w\ ■ = wx Ww^ = -w^Ww\ = -Aw^ • w\ , => (fi + \)wft ■w\ = 0. (C.3) 

The vectors and w\ are orthogonal (with respect to the scalar product w M ■ w\ = 
^l2j( w fj.)j( w ^)3j that is not the standard one because of the absence of complex 
conjugation) unless p = —A. If the eigenvalue A (and so —A) is rt-times degenerate, 
in general one can redefine the eigenvectors in the same n-dimensional eigenspace so 
that each vector becomes orthogonal to the other n — 1 ones. Thus the following 
decomposition holds 

W= Y A^^eV^!, (C.4) 
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where {A}/± means that each couple ±A contributes only once in the sum and the 
symbol v (£> w identifies the matrix with elements ViWj. By introducing the auxiliary 
matrices 

+ _ w\ (g) w-\ + w_a <8> w\ 



n 



W\ ■ W- X 



which satisfy together with II A the orthogonality relations 







Ha 












^+A n + , 



any function of W expresses in terms of 11^ 

f( w ) = E (/ e (A)n+ + / (A)n-) = 

{A}/± {A}/± 

with f e and f a the even and odd part of /. 



f(\)w\ ® U?_A + f(-X)w-x ® u;a 



u>a • W-A 



(C.5) 



(C.6) 



,(C7) 



The decomposition (C.4) allows to express in a factorized form 

-A [w x ■ a,w-x ■ a] 



Z(W)p w = exp( 



A [wx ■ a, W-x ■ a 



2 2«;a • W-x 
{A}/± A A {A}/± 



n h h ( 

{A}/± 



sinh 



11- 6XP ("2 2/r A • w_ x 
An [wx ■ a, w_a • a] 



2w x ■ w-x 



where we made use of w?a's orthogonality (|C.3|) 

= 0, 



Fa • a, w_a • aj , ^ • a, w_ M • 



u?a • a, w-x ■ a] = 4(w A • W-xY 



We obtain Z{W) by imposing the normalization of pw 



Z(W/) = Trexp(^) = JJ 2 cosh(^) = ±^det|e" + e " 



(C.8) 

(C.9) 
(CIO) 



{A}/± 



The correlation matrix Y%j = <5y — (ajOj) comes directly from the quadratic part of 
expression (C.8) 

/a T Wa\ 

X 



Tr— '- — exp( 



J 



ncosh(^) J2 tanh(|)Tr 

{A}/± {A}/± 

[n*«- ®Mt) 

{A}/± 



Qj)QiJ [wx ■ a, w-x ■ aj 
2 2w A • w-x 



that is 



r = tanh 



(?)■ 



(C.ll) 



(C.12) 



The correlation matrix does not determine univocally function Z(W). Indeed we 
obtain the expression 

2 n 



Z(W) 



det 1 - T 2 



(C.13) 
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with an evident phase ambiguity (that is actually a sign ambiguity). In spite of the 
ambiguity in Z(W) ~ Z[T], we believe that the trace of two fermionic RDMs 



Z{\og(e w e w ')) 



Tr[pwPw>] = 
is a functional of T 

{r,r'} = Trp[r] P [r' 



Z{W)Z{W) 

tanh(W/2) and F = tanh(W'/2). In fact, we argue that 
1 

2^ 



(C.14) 



(C.15) 



n a 

{A}/2 

where the product is extended over all eigenvalues A of 1+IT" with degeneracy (always 
even) reduced by half. We have not a rigorous proof of Eq. (C.15 1, but we give an 
argument to justify it. After infinitesimal variations of T and V , function Z(W) 
changes in the following way 



dlogZ(W) = ~Tr^^ 



thus from equation ( |C. 13 ) it follows 



j i m i T r x r'd(r x rp 

dlogTr pw Pw , = - Tr i - ^ 



TdT 



T'dV 



i -r 2 i 



(C.16) 



(C.17) 



Using Eq. (44) for T x T', many simplifications occur, leading to 

d(rr') 



:dr 



1 / 1 
rdr' = - Tr 

-IT' 2 1 



(C.18) 



dlogTr[ ^' ] = 2 Trr 'rTrr'^ ■ i + nv™ 2^1 + rr' 

The variations dr and dr' are skew-symmetric by construction. Our task is to 
integrate the differential equation above. We consider all possible sets of 2N smooth 
simple curves joining the IT' eigenvalues Ai from IT' = (when Tr pwPw — 2 _Ar ) to 
our desired values. Despite of the initial ambiguity in the sign of Z(W), the expression 
for the logarithm 

1 2N r dA 

] gTt\p wPwt ] = -J2 / I - j -l--.2V]og2. (C.19) 

i— 1 ^ i 

is well defined, in view of the fact that IT' has a double-degenerate spectrum and the 
ambiguity related to the choice of the curves 7, (avoiding any singular point) gives 
simply an additive constant proportional to 2iri. Considering the path characterized 
by d(rr') = dA(rr' — l)/2, the product rule can be expressed in the following integral 
form 



{r,r'} = ex P 



Tr 



IT' - 1 

2 - a + Arr' 



-dA 







f 1 


= exp 


[2 J 


1 + A 



rr'-i 

Tr dA 

1 + Arr' . 



7[o->i] 7[0-H] 

where 7 is a smooth simple path in C which joins to 1 and avoids any of the 
points at which 1 + Arr' fails to be invertible. Because of the factor 1/2 in front of 
the integral we get equation (C.15). For this proof we used several properties about 
Pfaffian reviewed in Ref. |38j . Summarizing 

P[TW'] = cxp(i J ^ Tr ^dA)p[r x T'] . 



(C.20) 
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A final remark about the XY correlation matrices in the double block case 
Ti = Fa and T 3 = Ta u b 1 /Tb 1 ( |Appcndix A ). The quantity 



(C.21) 



can be 1 or —1 and numerical data suggest that it is independent of the subsystem 
and of the chain length. In particular we find 

\h\ > 1 or 7 
otherwise, 



1 

-1 



0. 



(C.22) 



so that 



{fi, T 3 } - sign(|/i| - 1) + 2(J 7O 0(1 - \h\) 



det 



i + rir 3 



(C.23) 



We close this appendix with some speculations and some simplified (wrong) 
expressions to which one could arrive if treating too naively correlation matrices. 
If we would not worry about the signs and non-invertibility, we could write directly 
the result of a generic product like Eq. ( 48 ) as 

{ri,---,r„} 



det 



n: 



,Wi 



\ rniidet 



\ 



i 1 T \ 1 

nH^IHiii 



■r, 
r, 



(C.24) 



This expression is generally wrong. When substituting the chain correlation matrices, 
(1 — IY3) is usually non-invertible (Fig. CI). Thus, except for S2, equation (C.24) is 
wrong. 

By using this wrong expression one can try to see if any simplification occurs in 
the general Renyi entropy. In general it is not possible to push this approach, unless 
taking further wrong assumptions. Assuming for example the (wrong) hypothesis of 
commuting T matrices, after trivial algebra we have (the wrong expression) 



1 



log 



Cl-Ca »=1 



\ 



I 01 1 
detJJ- 



n 



i + r Ci 



i=l 



(C.25) 



We wrote this expression here, because one could have wonder whether it (despite 
wrong) would have been defined for any real a (i.e. depending only on the eigenvalues 
of the matrices T) to have an expression analytically continuable to a 1 (and maybe 
correct in this limit). However, even this simplified wrong expression is only calculable 
for integer a and thus we did not insist in this direction and we preferred to exploit 
the exact solution without assumptions. 



Appendix D. Accidental Singularities 

When the correlation matrix of the region interposed between the blocks Tb 1 b 1 is 
non-invertible, that is the o~ z string mean value vanishes (S) = 0, some problems 



arise. The density matrix ( 32 ) is not well defined and all successive discussions are 
meaningless. Here we give some tricks to solve these problems. We consider a generic 
operator 

oeA 
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Figure CI. Sorted absolute value of the eigenvalues of Ti and T3 for two identical 
blocks of length i\ = t% = 100 at the distance r = 50 of an infinite Ising chain. 
n is the eigenvalue's positioning number. The eigenvalues of Ti lie between — 1 
and 1 while the eigenvalues of T3 can be higher (lower) than 1 (—1). Almost all 
eigenvalues are in the neighborhood of 1 and —1. 



If (e^> £ 0, as we observed earlier, Qw is proportional to an effective free density 
matrix 

aWa 

Qw = (e — )p r[w] . (D.2) 

But when (e^r 5 ) vanishes (it can happen only when W is not hermitian) the previous 
expression is meaningless. The trick is to control the singularities introducing the 
traceless string operator Se = I\j<=E i^ij-i^j {S% = 1) that belongs to the subsystem 
Ed A 

Qw^S E ^Y. (e^O)S E =S E ~Y1 (^SeO)O = (e 1 ^ S e )Sep[T] , (D.3) 
oeA oeA 

and such that the mean value (e^T^ Se) is different from 0. The correlation matrix 

r = i- (e /if^ a) (D.4) 

(e 4 s E ) 

can be used to perform any product at the price to have the string Se- We observe 
that 

S E p[T] = p[P E TP E ]S E , (D.5) 

where Pe is the diagonal operator with elements i = — 1 if i € -E and 1 otherwise. 
In fact, exploiting the commutation relation, products of operators like Q\y can be 
carried to the form 

mm m 
i— 1 i—l i—1 
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where is the final a z string. The trace 

mm m 

Tr y[ = n< e ^^> Tr n<°[ r ^ = 

i—1 i—1 i—1 

m 

= H(e — S Ez ){r u ---,r m }Tvp[r 1 x---xT m ]S & (D.7) 
is free from pathologies. 
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